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1. Introduction. 

In recent years it has become increasingly clear that a problem of central importance for 
the understanding of disordered spin systems is the control of random fluctuations of ther- 
modynamic quantities [AW,NS,BM,T1]. Unfortunately, a precise control of such quantities 
is very hard to come by. Concentration of measure techniques [T2] have been realized to 
be efficient tools to get upper bounds [BGP1,BG1], but lower bounds or exact limit theo- 
rems are scarce. One of these examples is the Sherrington-Kirkpatrick (SK) model in the 
high-temperature phase, where a central limit theorem for the free energy was proven first 
by Aizenman, Lebowitz and Ruelle [ALR], using cluster expansion techniques, and later by 
Comets and Neveu [CN], making use of martingale methods and stochastic calculus. Their 
methods have been extended to a few related cases [Tou,Bl] later. In the present paper we 
want to continue this effort by investigating a large class of natural generalisation of the SK 
model, the so called p-spin SK models, and their p | -|-oo limit, the random energy model 
(REM). 

For our present purposes it is natural to consider the class of models we study as Gaussian 
processes on the hypercube Sn = {—1, 1}^- We will always denote the corner of Sn by a; 
for historical reasons they are called spin configurations. A Gaussian process X on Sn is 
characterized completely by its mean and covariance function. The processes we consider 
will always be assumed to have mean zero and covariance 

EX^X^, = f{RN{a,a')), (1.1) 

where / depends on the so-called overlap^, Rn{(t,o'') = N~^{a,a') = N~^J2iLi(^i(^'ij 
this note we will concentrate on the case where f{x) = fp{x) := x^, with p even.^ In this 
case, Xa can be represented in the form 

11,12, ■■■ ,ip 

with Ji^,...^i^ a family of i.i.d. normal random variables. Since iorp = 2 we obtain the classical 
SK model, this representation justifies the name p-spin SK model. Note that as p increases, 
the process gets more and more de-correlated, and in the limit p | -|-oo we arrive at the case 
where X^- are i.i.d. normal random variables. 

^The overlap is related to the Hamming distance dnam by dHamicjc') = -^(1 — J?jv(f, <t'))/2- 
^The case p odd can also be treated, but presents considerable additional computational problems. 
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Given such a Gaussian process, our main object of interest is the so called partition func- 
tion, 

Zis,N ^ E.e^^^" ^ 2-^ e^^""^- (1-3) 

The quantities qP^^'^ arc called Boltzmann weights and the parameter j3 € M+ is known 
as the inverse temperature, and Hf^{a) = VNX„ as (minus) the Hamiltonian in statistical 
mechanics. ^/3,7v are random variables, and we arc primarily interested in their behaviour as 
tends to infinity. In statistical mechanics, it is customary to introduce the so-called free 
energy 

Fp,N = -^lnZ0,j,. (1.4) 

It is easy to prove in all the models we consider here, that for all values of P, F^^n is a 
self-averaging quantity, i.e. that 

lim |F^,Ar -Ei^)3,iv| = a.s. (1.5) 

It is, however, not known in general whether the so called quenched free energy EF^^^r con- 
verges to a limit as N tends to infinity. This has, however, been proven for sufficiently small 
values of /3: more precisely, one knows that 

Theorem 1.1. Define ^2 = '^, and for p > 2 

inf (l + m-^).^(m) (1.6) 

0<m<l 

where <p{m) = [(1 — m) ln(l — m) -|- (1 -|- m) ln(l -|- m)]/2. Then for all (3 < Pp 

lim = -/?/2. (1.7) 

Af 1+00 

Remark: For p = 2 this result was first proven in [ALR]. A very simple proof has later 
been given by Talagrand [T]. Comets [C] has shown that the value /3 = 1 is optimal in the 
sense that (1.7) fails for /? > 1. The result for p > 3 is due to Talagrand [Tl]. It is clear 
that in all cases (1.7) will fail for (3 > V21n2 which by a more elaborate computation can 
be improved to /3 > \/21n2(l — 2~'^'p'p) with Cp < 5, for p large [B2]. On the other hand, a 
simple calculation shows that ~ \/2 ln2(l — 2~^'/21n2). One should note that to get (1.7) 
up to a value so close to \/21n2 required a substantial modification of the original argument 
of [T2], namely the use of a "truncated" second moment method. Such a truncation will also 
be the main difficulty in obtaining our results^. 



^For similar reasons, slightly diflferent truncations were also used by Toubol [Tou] (and probably first) in 
the study of the CLT for the SK model with vector valued spins. 
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In the case of the REM, it is well known that the critical inverse temperature (3rem 



V21n2 and that [D2] 

/ -/3/2, if/3<V2T^ 

lim r n Af REM — S / , / 

iVT+oo I -V2hr2 + /3-Mn2, if /? > ^21^. 



As a consequence one has that (this result is essentially contained in [Dl], a rigorous proof 
follows easily from the results contained in [Tl]^) 

lim limsupF/3,Ar,p = lim lim inf F^,Ar,p = lim Fp^N,REM- (1-9) 

In this note we will control the fluctuations of the free energy in (essentially) all of the domain 
of parameters /3, p (even) where the limit is known to exists, i.e. the high temperature regions 
of the p-spin models, and the entire temperature range in the REM. Although the REM is 
rather singular and the techniques used for that case are totally different from those we will 
use for the p-spin models, we felt it would be instructive to include this singular limiting 
case in this paper. Moreover, it turns out that in spite of the heavy investigation the REM 
has enjoyed over the years [Dl,D2,0P,GMP,Ru], no precise fluctuation results for the free 
energy are available in the literature. Finally, we are convinced that the reader will be rather 
surprised by the rich structure the fluctuation behaviour this model exhibits. 

Let us now state our results. We begin with the j3-spin SK models. 

Theorem 1.2 Consider the p-spin SK-model with p = 2k > 2. There exists (3p > such 
that for all P < (3p 

in distribution as N 1 +oo, where Moo(t) is the centered Gaussian process with mean zero 
and covariance 

E{M^{t) - MM? = {t-s) [{p - 1)!!] . (1.11) 

The value of jSp can be estimated reasonably well. To state lower bound on f3p we need, 
however, some notation. We define the functions 

/(mi, 7712, ms) = ^((1 + mi +7712 + m3)ln(l +mi +m2 + ms) 
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+ (1 — mi — m2 + ms) ln(l — mi — m2 + mz) 
+ (1 + mi — m2 — ma) ln(l + mi — m2 — ms) 
+ (1 - mi + m2 - ms) ln(l - mi + m2 - ms)^ , 



(1.12) 



^Private communication by M. Talagrand. 
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5p(mi,m2,m3)= (l + ^ + 1 + ^ ^ 

L ^ 1 — TTlg ' ^ 1 — 



P P > 9 

1- Tr7 „ \ 



+ 



2mf( 



1P 



p^p 



1 + 



1 -mg^ 



Km5 — m%m%, \ 
1 + ^ — If^J 
1 — ' 



1/2 



np{mi,m2,'m3) — — 



2p 



2(1 - mg 



2pN 



(1.13) 



(1.14) 



and 



(1.15) 



Up{mi,m2,ms) =I{mi,m2,ms){l + mg) S'p(mi, m2, 777,3)^2 + 2777,3 

+ i?p(mi,7TT.2,77i3)(l + mf ) — (2 + ) 

on the set 

A ={?7ii, m2, ms G [—1, 1]^ I 1 — mi — m2 + m3 > 0, 1 — mi + m2 — m3 > 0, 

1 + mi — m2 — m3 > O}. 

Note that the function 7(mi, m2, ms) is symmetric in mi, m2 and ms, and that ^^(mi, m2, ms), 
Rp{mi,m2,m3) and Up{mi,m2,mz) are symmetric in mi and m2. Let 



;i.i6) 



I^(mi,m2,m3) =max-^ /(mi,m2,m3) ( - + 



3 ' m 1 + m2 + mg y ' 

Up{mi, 7712, 1713), Up{7ni, 7713,7712), Up{77l2,77l3,77li)y 

With this notation we have 
Theorem 1.3 Let p = 2k > 2. Then 



(1.17) 



inf yp(mi,m2,m3) < /3p < 21n2. 

mi,m2,m3EA 



In particular 



lim (5; = 21n2. 



(1.18) 



(1.19) 



We see that the scale on which the partition functions fluctuate decreases rapidly as p 
increases. One might guess that the scale becomes exponentially small in N in the limiting 
random energy model. This is indeed true, but more surprising things happen, as the following 
theorem states: 

Theorem 1.4 The free energy of the REM has the following fluctuations: 
(i) If/3< ^JWlj2, then 



(1.20) 
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(ii) If 13 = A/In 2/2, then 



^^(0,1). (1.21) 



(in) Let a = /3/V21n2. // i/ln2/2 < /3 < V21n2, i/ten 

-(^/2h^2-/3)2 + f [ln(Arin2)+ln47r] ^/3,iV 



/oo 
e«^(7'(dz) -e-^dz), (1.22) 
-oo 

where V denotes the Poisson point process on TZ with intensity measure e~^dx. 



e 2 

EZ^,jv 



Theorem 1.3 covers the high temperature regime. However, in the REM we can also 
compute the fluctuations in the low temperature phase. 

Theorem 1.5 

(i) If 13 = V21n2, then 

^1 [in(iv in 2)+in 4.] f ^ _ 1 ln(iV In 2) + In 47r X ^ j T.p^^,). 



(1.23) 



(ii) If f3> V21n2, i/ien 

oo 

g-iV[/3V2hr2-l„2] + f [l„(iVl„2)+ln4.]^^_^ ^ y (1.24) 



and 



In 



OO OO 

Z^,iv -ElnZ^,iv ^ In j e"^7'(dz) - Eln ^ e'^^V{dz). (1.25) 



Remark: Note that expressions like f^^e""^ {(P{dz) — e~^dz) are always understood as 
liuiyi-ac Jy G°"^ {'P{dz) — e~^dz). We will sec that all the functionals of the Poisson point 
process appearing are almost surely finite random variables. 

Remark: Note that the Poisson integral e"'^V{dz) is the partition function of Ruelle's 
version of the REM [Ru]. Thus (1.25) affirms that above the critical temperature, the fluc- 
tuations of the free energy of the REM converge in distribution to those of Ruelle's model. 
While this connection was surely evident for Ruelle and motivated the introduction of his 
model, we have not been able to find a rigorous statement of this connection in the literature. 
In [GMP] the scale on which fluctuations take place has been established, but no actual limit 
theorem was proven. 
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Remark: It is interesting to observe that in the REM there is a second "phase transition" 
within the high-temperature phase at which the fluctuations become non-Gaussian. In fact, 
in the REM the main phase transition can be interpreted as a breakdown of the Law of 
Large Numbers, while the second transition corresponds to a breakdown of the Central Limit 
Theorem. 

The remainder of this paper is organized as follows. In the next section we present the 
proofs of Theorems 1.2 and 1.3. They are based on an adaptation of the martingale method of 
Comets and Neveu. The essential new ingredient is the rather involved truncation procedure 
inspired by Talagrand's work. However, in the proof of the CLT, the computational aspects 
become even more involved and require the consideration of truncated third moment of the 
partition function. For this reason Section 2 is rather long and quite technical. However, the 
proof is organized in such a way that the CLT is first proven for "very high" temperatures 
where no truncations are necessary, while the more technical aspects needed to approach 
the critical temperature are dealt with separately later. Section 3 is devoted to proving 
Theorems 1.4 and 1.5 for the REM. It is technically completely different and independent 
from Section 2. It can therefore be read independently from the rest of the paper. In an 
appendix we explain some of the technical difficulties that appear in the case p odd and we 
explain the result to be expected in that case. 

2. The CLT in the p-spin model 

The proof of the central limit theorem in the p-spin SK model relies on a martingale 
central limit theorem which uses that fact that a Gaussian random variable can always 
be seen as the marginal distribution of a Brownian motion. Thus wc follow Comets and 
Neveu and introduce the p-parameter family of independent standard Brownian motions 



iJiui2,...,iA^),t G K+)n..2,...,*peN withEJ,,,i„...,i^(t) = and E ^^^^^^^ (t) = t. The Hamil- 



Note that we can also consider it as a Gaussian process on {-1,1}^ X M+ with mean zero 
and correlation function 



where fp{x) = x^. In particular, we have KHp^la, t) = Nt and E eyi]i{Hi\]{t, cj)} = exp{A''i/2} 
for all a. For later convenience we introduce the normalized partition function 




(2.1) 




(2.2) 



ZM{t) = E^ exp{i/jv(i, a) - Nt/2}, 



(2.3) 
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It is related to the partition function ^/3,jv of Section 1 by Zn{P'^) = Z^^jv/EZ^.a?, with 
equality holding in law. The important point of this construction is that that for all fixed 
N > 1, Z^it) is a continuous martingale in the variable t with EZjv(t) = 1. 

We begin the proof with some preliminary steps along the lines of [CN]. Let us find the 
bracket < Z]s[{t) > of the martingale Z^it), i. e. the unique increasing process vanishing at 
zero, such that Zjf{t)— < Z^it) > is the continuous martingale (see [RY]). By Ito's formula, 
ZN{t) satisfies the following stochastic differential equation: 

dZiv(i) = exp{ifiv(i, a) - Nt/2} di?iv(i, o"). (2.4) 

Then due to well-known properties of martingale brackets 

t t 
< ZN{t) >=E^,a' < Je^"^''"^-^'^^ dHN{s,a), Je"''^''''"^-^'/^ dHN{s,a') > 



t 

=E^,^/ J e-f^^(«.<^)+-f^^(«'^')-^« d < Hn{s, a),HN{s, a') > (2.5) 



t 

=E,,,. J e^-('*'")+^-(^'"')-^''iV/p(i?jv(c7,a')) ds. 



Since 

E / Z-^(s) d < Zn(s) >=E / „ , , \, ^ ds<Nt< oo, 

^'^ 

(2.6) 

we may introduce a continuous local martingale MN{t) = Z^^{s) dZN{s). Thus Z^it) 
solves the stochastic differential equation 

dZuit) = ZN{t)dMN{t) 

and the following fundamental representation of Z^it) holds: 

ZN{t) = exp{MAr(t) - 1/2 < Mjv(t) >}. (2.7) 

Here < M^it) > is the bracket of MM{t) and < MM{t) >= Z^^{s) d< Zm{s)> . Let us 
note that 

^<M^{t)> = Z^\t)^<ZN{t)> 

at at ^2.8) 



= Z/(t) (^E.^^,e^-(*'-)+^-(*'-')-^* iV/p(i?iv(a, a'))) • 
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Note also that Mj>f{t) is locally square integrable. In fact, by (2.6) 

t 

EM^(t) = E < Mjv(t) >=^J ^jv^(«) d < Zn{s) >< Nt < oo. (2.9) 



To prove Theorems 1.2 and 1.3, we will show that for all t satisfying 

t< inf Yp{mi,m2,'m3). (2-10) 

the bracket of the local martingale Ar(p-2)/4Mjv(i), which is Ar(p-2)/2 < Miv(t) >, converges 
to tE^P in probability as AT f +oo. Here ^ is a Gaussian random variable with = 0, 
E^^ = 1. Then by the martingale convergence theorem (see Theorem 3.1.8 in [JS]) the local 
martingale Ar(P~2)/*Mjv(t) converges to M^oit) in law as AT | +oo. This fact together with 
the representation (2.7) implies immediately the statement of Theorem 1.2. 

Sketch of the proof of Theorems 1.2 and 1.3: We will now outline further steps of 
the proof. First, we show the convergence on a more restricted 

interval of t. Lemma 2.1 reduces this problem to the convergence of 

Ar(P-2)/2E|1/^(i)| ^ 0, asiVt+oo, (2.11) 

where 

VN{t) := Ar-^E^,<,.(7VP/2/p(i?iv((T,(T')) -ECP)e-f^'^(*''^)+^^(*''^')-^*. 
The proof of this lemma is based on the fact that 

^(p-2)/2 d ^ ^^^^^ ^ _^^p ^ ^(p-2)/2 M^l^ (2.12) 

at Zpfit) 

and is performed via integration. It almost mimics the proof proposed in [CN]. In particular, 
we use the fact that Zlf{t) is not small on events of large probability. The convergence (2.11) 
is proved in Proposition 2.2. Let us give some intuition for it. One can write 

EVNit) = (Nfpim) - iV(2-p)/2E^P)e*^^''("^)p(c7 • a' = mN). (2.13) 

m=0,±l/Ar,...,±l 

By Stirling's formula 

P((7 • u' = mN) ^ e-^-^^"^), 

^J2^^{1 + m){l - m)N 
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where (f){m) = [(1 + m) ln(l + m) + (1 — m) ln(l — m)]/2. (here and in the sequel we use the 
symbol ~ to denote asymptotic equivalence, i.e. qn ~ ^at limAr^+oo = !)• Note that 
(/>(m) = — m^/2(l + o(l)) as m ^ 0. Now split the right-hand side of (2.13) into two terms: 
the summation in the first one will be over m with \m] "small enough" and in the second — 
over all other m. It is not difficult to treat the first term. Since p > 3, we have for any fixed 
t 

tfpim)+(j){m) = -m'^/2{l+o{l)), m ^ 0. (2.14) 



Then putting my N = s, the first term becomes 

oo 

9 7\r(2— p)/2 f 

-^L= V (iv(2-p)/2gP _ jv(^-P)/^E$P)e-^'/" ~ ^ / (sf - Ee)e-^'/2 ds, 

s=m\/N —oo 

from where the normalisation Ar(P~^)/^ is immediate. To ensure the convergence to zero of 
the second term (the one with correlations m not close to zero), the power of the exponent 
in it should be negative: 

sup {tfp{m) - 4>{m)) < 0. 

me [0,1] 

Thus for all t < info<™<i (t){m)m-P , we get Ar(p-2)/2EV;v(t) ^ 0. Note that. Proposition 2.2 
states a stronger result (2.11). To get rid of the absolute value of Vjv(i) in (2.11), we follow an 
idea suggested in [CN] to apply the Cauchy-Schwartz inequality. Thus, instead of E|y;v(i)|, 
we get WN{t) (see the proof of Proposition 2.2) which refers to the third moment of ZN{t). 
This makes technical computations slightly tougher and leads to the bound on t (2.19) given 
in Lemma 2.1 below. 

Note also that these arguments arc valid only for p > 3. The case p = 2 of [CN] and 
[Tou] is different, since there, (2.14) does not hold. This case is treated in [CN] by the 
multi-dimensional Central Limit Theorem for N independent vectors (cTjO"^, a^cr", (Tjcr"). 

Next, we will extend the bound (2.19) to the full regime announced in (2.10). We have 
seen, that (2.19) was imposed by configurations of spins with rather big correlations m in 
the sum (2.13). We will reduce their contribution, using Talagrand's idea to truncate the 
Hamiltonian. Consider instead of VN{t) 

VN{t,e) =E^,^,(^Nfp(^RN{a,a')^ _ jy{2-p)/2^^pyH^(t,a)+H^{t,a')-Nt 

X '^{HN(t,a)<(l+e)tN,HNit,CT')<{l+e)tN} 
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for some e > 0. Then 



EVn {t,e)= J2 i^fp ("^) ^ tN^^~^^/^ )F{a-a' = mN) 

m=o,±i/Ar,...,±i (2.15) 



^ -^{^l<^/JVt(l+e),6<^/JVt(l+e)}' 

where ^1,^2 are standard Gaussians with cov(^i,^2) = fn. Let us again split W/j^{t,e) into 
two terms with "small" and "large" m in the sum (2.15). The analysis of the first term is 
completely analogous to the one in the case of VN{t). Wc can neglect the truncation here, 
since ^1 and ^2 are almost independent. In the second term, ^1 and ^2 arc more correlated. 
But due to the truncation, the expectation of the exponent involved in this term is much 
smaller than e^'^''^ . In fact, by the elementary estimate (5.2) for Gaussian random variables 

117 v/iVt^l+ViVt^2-JVt-,r 

-^{€i<V^(l+e),S2<V^(l + e)} 



Then for any 



< IP yArt(2+2mf)^-Artir 

- -^{^<2v^(l + e)(2+2mP)-i} 

< exp{[-4Ari(l + e)2][4 + 4mP]-^ + 2Nt{l + e) - Nt} 
= exp{[ArimP(l + 2e) - Nte\l + mP]-^}. 



t< inf {l + m-P)(p{m) (2.16) 

0<m<l 



and for an appropriate choice of e all terms of the sum (2.15) with m not close to zero are 
exponentially small. This implies N^P~'^'>/'^'EVN{t, e) 0. The bound (2.16) is Talagrand's 
bound for the critical temperature in the p-spin SK model, see (1.6). It tends to 2 In 2 as 
p T +00. 

In order to incorporate this idea into our proof, we reduce the problem of convergence 
jY(p-2)/2 ^ MN{t) >^ iE^P to the following statements: 

7V(P-2)/2E|y;v(«,e)| ^0, (2.17) 

and 

iVb-2)/2E|(VV(t) - VN{t,e))Z-\t)\ ^ 0, (2.18) 

for all e > 0. This is derived in Lemma 2.3 again from (2.11). In Proposition 2.4 we show 
(2.17). Again, because of the absolute value, we must apply the Cauchy-Schwartz inequality 
and pass to the third moment of Zj^it). This makes technical computations much harder. 
Namely, we get three standard Gaussian random variables ^1,^2,^3 with covariances mi, m2. 
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ma- To benefit from the truncation for obtaining a good bound on t, wc have to take into 
account four different cases: one when all mi, 1112,1713 are large and others when two of these 
correlations are large and the third is small. Then the analogue of (2.16) is the minimum of 
four estimates of this kind. Therefore, the bound (2.10) is the minimum of four functions. 
The convergence (2.18) is the subject of Proposition 2.5. Its proof uses ideas of M. Talagrand 
and a concentration of measure inequality. 

Lemma 2.1: Let 

T<mf ^<;";- ";'■""> . (2.19) 



Then 



sup |Ar(p-2)/2 < Miv(t) > (2.20) 

0<t<T 



in probability, where ^ is a Gaussian random variable with = 0, E^^ = i 
Proof. Let us denote by 

VN{t) = ^ < z^it) > -N(^-py^z%{t)Ee. 

Then 

^^(p-2)/2 ^ ^^^^^ ^ _^^p =N^P-^y^VN{t)Z^^{t) 

=N(P-^y^VNit)exp{-2MN{t)+ < MN{t) >}. 

Let us introduce the events 

A^t^ ■= {-Mjv(0 < a + (6/2) < Mjv(t) > for all t > 0}. 

Note that by an appropriate choice of a > and 6 > 0, their probabilities can be made 
arbitrarily close to 1. In fact, the process -Bjv(t) = Mjv(5't), where St = min{s |< Mn{s) >= 
t}, is a standard Brownian motion and M^it) = Bn{< M^it) >). By the well-known fact 
for Brownian motion 

P{^^J = P{-Siv(t) < a + {h/2)t for alH > 0} > 1 - exp{-a6}. (2.21) 

We have: 



= iV(^-')/'|Fiv(t)|exp{-2Mjv(t)+ < MM{t) >}W } (^.22) 
<jY(p-2)/2exp{2a}|yAr(t)|exp{(l + 6) < Mjv(i) >}. 
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Let US also introduce the function Xb{x) '■= [1 — exp{(l + 6)a;}][l + b] ^. Then by (2.22) for 
alH < T 

|7V(p-2)/2^^(< Mjv(t) > -iAr(2-p)/2E^P)]i I 

t 



X exp{-(l +6)(< Mjv(s) > -siV(2-p)/2E^P)}ds 

t 

<jY(p-2)/2 /■ < Mjv(s) > -iV(2-P)/2E^P I.^jv . 

7 ds 



X exp{-(l + 6)(< Mjv(s) > -sAr(2-p)/2E^P)}ds (2.23) 

t 

<Ar(P-2)/2exp{2a + TiV(2-p)/2(i + 6)} ^ |V;v(s)|ds. 



This yields 



iV(P-2)/2 sup |x6(<M^(t)>-iAr(2-p)/2E^P)|i ^^ 

(\^ j.^ ri-f ^ a. 



0<t<T 



(2.24) 



<Ar(p-2)/2exp{2a + riV(2-p)/2(i + 5)} ^ |y^(s)|ds. 



We will show in Proposition 2.2 that 



lim N^P-^^/^E\VN{t)\ = 

N1+00 

uniformly in t € [0,T]. Consequently supjv>i^t<T -^'■^~^'*^^lE|Vjv(OI < Then by the domi- 
nated convergence theorem 



lim E\N^^ ^•'/^ sup \Xb{< MM{t) > -tiV(2-p)/2E^f)|I. 



0<t<T 



It follows that for all a, 6 > 

^(p-2)/2 g^^p 1^^^^ ^^^^^ -tAr(2-p)/2E^P)|I.^jv . ^ as AT t +00 



0<t<T 



in probability. Then also N^^ 2)/2 supg^^^^^ 

|X6(< Mjv(t) > -tiV(2-p)/2E4P)| ^ oo, as 
iV — >■ oo since by (2.21) the probability of the events A^^^ can be made arbitrarily close to 1. 
This last fact imphes (2.20) and the lemma is proved.<0> 



It remains to prove the following proposition. 
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Proposition 2.2:Assume that T satisfies (2.19). Then 

lim Ar(P-2)/2E|VAr(t)| =0 (2.25) 

iVt+oo 

uniformly in [0,r]. 

Proof. It follows from (2.5) and the definition of Fjv(t) that 
By the Cauchy-Schwartz inequality: 

7V(P-2)/2E|^^(t)| =E|E,e^^(*''^)-^*/'E,, (n^'^ u{RN{a,a')) - ^^pyH^{t,a')-m/2 
< EE,e^^(*''^)-^*/2l'^' 

'^a,a',a" (iV^ / (i?iv (a, a')) - E^) (iV^/ (^AT (a, a")) - E^) 

X exp{iVt(/p(i?Ar(a,(7')) +/p(i?iv(<T,0) +/p(^iv((T((7',(7")))}" 
Then it suffices to prove that 

X exp{Art(/p(i?^((7,(7')) +/p(i?iv(cT,a")) +/p(i?iv(a(a',a")))} 
tends to zero uniformly in [0,T] as iV | +oo. We represent it as 

mi,m2,m3e.Ajv 

X P{(7 • o-' = miN, u-a" = m2N, a' ■ a" = nisN} 

where the set An = ACi {0, ±1/N, ±2/N, . . . , ±1}'^. A standard combinatorial calculation 
yields 

F{a ■ a' = miN, a ■ a" = m2N, a' ■ a" = msN} 



1/2 



= 2 



-2N 



N{l-mi)/2 



(2.26) 



N W A^(l + mi)/2 
yN{l + mi)/2y \N{1 + mi +1712 + m3)/4:J ViV(l + ma - mi - m3)/4) ' 

By Stirling's formula we obtain 
P{cr . a' = rriiN, a ■ a" = maiV, a' ■ a" = m^N] 

16exp{-Af/(mi,m2,m3)} , , \n _l m-i/2 

= ^(27r)3A^3 [(l + mi + m2+m3)(l-mi-m2 + m3)] ' ^2.27) 

X [(1 + mi - m2 - m2)(l - mi + m2 - ms)]"-^/^^! + oi^-^^ as N ] +oo, 



Fluctuations 15 

for any given m-\_,m2,ms E An- Let us remark that 

t(m? + m| + ml) + (m? +ml + ml)/2 - 7(mi,m2,m3) = 0((|mi| + |m2| + Imgl)^) (2.28) 

as mi, m2, 1713 uniformly in [0, T]. Then for all sufficiently small e > there exists a 
constant h> such that 

sup [t{ml + m2 + m^) — I{mi,m2,m3)~\ < —h{ml + m2 + m^)/2 (2.29) 
te[o,T] 

for all mi,m2,m3 e ACi {\mi\ + \m2\ + |m3| < e}. Let us fix such a small e > and an 
arbitrary constant < (5 < 1/6 and then split Wjv(i) into four terms: 

where 



mi,m2 ,m3 EAj^ 
l">ll + l'"2l + l"'3l<-^~^^^~^ 



16 ^r^•™2 , „2 , /o\ 



jV(m?+m^+m2)/2 



/3^(t) = (^NP/^fp{mi) - E^) (iV^/Vp(rn2) - E<e^')e^*(-'^f 



1^1 1 + 1^2 1 + 1^3 l>-^~^/^~'^ 
l^^l H-|'"T-2 l + l'^2 l<€ 



X F{a ■ a' = miN, a ■ a" = m2N, a' ■ a" = m^N), 



l^ll + l'^2l + l^3l>^ 



X P((7 • a' = miN, a-a" = m2N, a' ■ a" = m^N). 
We will prove that all four terms l]^,l1f{t),I^{t),lff{t) tend to zero uniformly in [0, T] as 
N t +00. 

To show this for let us put mi^/N = si,m2\/iV = S2-im3\/N = S3. Then 

hm Ij,= hm ^ y (sP-Ee)(4-EOe"^''+'^+'^^/' 

= 0,±1/ v^,±2/\/jV,... 
|sil + |s2l + |s3l<-^^'''*"* 



OC' OC' 00 



^== I I y"(xP-En(yP-EOe-("'+^'+^')/Wydz = 0. 



-00 —00 —00 
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To treat I%{t), we rewrite it using (2.27) as 



\mi\ + \m2\ + \m3\<N-'^/^-l 

pA''[t(/p(mi)++/p(m2)+/p(m.3))+(mi+m|+m3)/2-/(mi,m2,m3)] 

X [(1 + mi + 777,2 + "i3)(l - mi - m2 + ms)]"^/^ 
X [(l + mi-m2-m3)(l-mi+m2-m3)]-^/2j|i + Oj^l.^^ _1 

(2.30) 

Moreover, here 0(1) is bounded uniformly in An H {|mi| + |m2| + |m3| < e} by Stirling's 
formula. It follows from (2.28) that 



lim sup iV|f(mf + mg + mf) + (mi + m2 + m3)/2 — 7(mi, m2,m3)| = 0. 

JVt+oo te[o,T] 

l'"ll + l"'2l + |m3|<Ar-l/3-'5 



Then 



lim sup 

N^+oo te[o,T] 

l'ni| + |m2l + |m3|<Ar-l/3-'5 



g-^[t(/p (™'2)+/p(m3))+('Tii+m2+m3)/2-7(mi,m2, 7713)] 



X [(1 +mi + m2 + m3)(l — mi — m2 + m3)] 

X [(1 +mi - m2 - m2)(l - mi + m2 - ms)]"^/^^^! + o(^-^ 



0, 



while 



lim 

JVT+OO , 



-77^ E I (^"^Vp(mi) - E^) (iV-/Vp(rn2) - E^^) 



-N{ml+ml+ml)/2 



lim 



l"lll + l"12l + l"13l<'V-^/3-* 

16 



^(27rA^)3 



^ \{sl- W) [si - W) 



«1 .»2 .■'3="'±l/v'^.--- 
|sil + N2l + N3l<JVl/«-* 



00 CxD CxD 
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J J J \{xP-W){yP-W)\e-^''^+y^+''^/^dxdydz<oo. 



-00 —00 —00 



Thus /^(t) ^ uniformly in [0, T] as iV j +00. 
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To estimate lff{t), we rewrite it in the same way using (2.27): 

\m-^\ + \m2\ + \m3\>N-'^/3-S 
\'mi\ + \m2\ + \mg\<e 

gAr[t(/p(mi)+/p(m2)+/p(m3))-/(mi,m2,m3)] 

X [(1 + mi + m2 + m3)(l — mi — m2 + ms)]"-*-/^ 

X [(1 + mi - m2 - m2)(l - mi + m2 - ms)]"^/^^! + o(^^ 

(2.31) 

Due to (2.29), there exists a constant h' > such that for all sufficiently large A'" 

sup exp{— Ar[i(mf + m2 + m3) — /(mi, m2, ms)]} 

te[o,T] 

\mi\ + \m2\ + \m3\>N-^/^-^ 

< sup ex.p{-Nh{ml + m^ + mi)/2} < ex.p{-h' N^/^''^^}. 

\m^\ + \m2\ + \m3\>N-'^/3-S 
\mi\ + \m2\ + \m3\<e 

The sum 

l'"ll + l'"2l + l™3l>-^^^''^~* 
Imi | + |m2l + |m3|<e 

has polynomial growth as A?^ t +oo and the uniform convergence I%{t) — in [0, T] is proved. 

Finally, let us consider I%{t). By Stirling's formula there exists a constant C such that 

for all (mi,m2,m3) G An n {|mi| + |m2| + \m^\ > e} 



P{(7 • a' = miN,a- a" = m^N.a' ■ a" = m3N}<CVNexp{-NI{mi,m2,m3)}. (2.32) 

Then by the assumption (2.19), for given T there exists a constant h" > such that 

sup exp{iVt(mf + m^ + mf )}P{c7 • a' = rriiN, a ■ a" = m2N, a' ■ s" = msN} 

te[o,T] 

\mi\ + lm2i + [mg\>e 



<CVN sup exp{-Ar[t(mf + m| + mf ) - /(mi, m2, mg)]} < CVNex.p{-h"N}. 

te[o,T] 

\mi\ + \m2\ + \m3\>e 

The remaining sum in this term has again polynomial growth, whence I%{t) — ^ uniformly 
in [0,r]. The lemma is proved.<0' 

Remark. Let us note that the restriction (2.19) on T was essential only for the analysis 
of the fourth term I%{t). This means that the convergence Ar(p-2)/2E|T/jv(t)| ^ breaks 
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down for larger T only because of the configurations of spins with rather big correlations 
a ■ a' = mi, a ■ a" = m2, cr' ■ <y" = ra^. To extend our result to the whole interval (2.10) 
of admissible T, we need to reduce the contribution of these configurations into Wjv(i)- For 
that purpose we will follow the idea of M. Talagrand [T] to truncate the Hamiltonian. 

Now we prove the statement of the previous lemma for all T satisfying (2.10). 

Lemma. 2.3: Let 

T< inf y(mi,m2,m3). 



Then 



sup I < Ar(p-2)/4Mjv(i) > -tE^P\ (2.33) 

0<t<T 



in probability. 

Proof. Let us fix e > such that for some constants /ii , /i2 > 

sup [t{mi + m2 + m^) — I{mi,m2,m3)] < —hi{mi + m2 + rri^) (2.34) 

*e[o,T] 

and 

sup t[mm{Qp(mi,m2,m3,e), Lp(mi, m2, ma, e), 



tG[0,T] 

7TX ^ + TTl + TTl ^ > 3 e 



(2.35) 



Lp{mi, 1713, 1712, e), Lp{m2,m3,mi,e)}] - /(mi, m2, ma) < -/12 

where 

Qp(mi, m2, mg, e) =[-9e'^ + 6(1 + 2e)(m? + m^ + mf )][2(3 + 2m? + 2m^ + 2mf )]-\ 



Lp{mi,m2,m3,e) = - 1 - mg - (1 + e)^ + (1 + e)Sp{mi,m2,m3)y 2 + 2ml 

+ i?p(mi,m2,m3)(l + m|) [l + mf]~^. 

Condition (2.34) is the same as (2.29) and, due to (2.28), for any given T > it is possible to 
find an appropriate e > such that (2.35) is satisfied. However, e > ensuring (2.35) does 
exist, if and only if T satisfies the assumption (2.10). The meaning of (2.35) will become 
clear in the proof of a further Proposition 2.4. Let us introduce 

VN{t,e) =E,,,.(Af/p(i?^(c7,a')) - A^('-^)/'EeP)e-f^^(*'")+-f^^(*'"')-^* 

X ^{HN{t,cr)<{l + e)tN,HNit,<7')<{l + e)tN} 

VN{t,e) =E,,,.(A^/p(i?^(a,a')) - A^('-^)/'Eep)e-^^(*'")+-^'^(*'"')-^* 

X I{l/jv(t,cr)>(H-e)tAr, Or HN{t,CT')>{l+e)tN} 
=VN{t)-VN{t,e). 



Fluctuations 



19 



Let us also fix some Tq > satisfying the assumption (2.19) of the previous lemma. Proceed- 
ing along the lines of the proof of Lemma 2.1, we get for all t e [To,T]: 



To 



< Ar(?'-2)/2 exp{2a + ToN^''-p^/^{1 + b)} J |Viv(s)| ds 



t 

+ N^P-^^/^e^p{2a + tN^''-P^/\l + b)} J \VN{s,e)\ds 

To 

t 

+ N^P-^y^ J \Vn{s, e)\Z^\s) exp{-(l + b){< Mj,{s) > -sN^^-p^/^)}!^^.^ ds. 



Then 



7V(P-2)/2 sup ^F^^^ MN{t) > -tiV(2-P)/2E^P)|]I.^jv . 
To<t<T 

To 

< iv(f-2)/2 exp{2a + roA^(2-p)/2(j ^ ^ |^^(^-)| 



T 

Ar(p-2)/2exp{2a + rAr(2-p)/2(i + t)} y" \Vj^{s,e)\ds 

To 

T 

iV(p-2)/2exp{riV(2-p)/2(i + 6)} y \VN{s,e)\Z]:^\s)ds. 



+ ■ 

To 

T 

^ o\ /o , f n ^\ /o , . - ^ ^ 

To 

It was proved in Lemma 2.1 that iV(P-2)/2E|VAr(t)| uniformly in [0, Tq] as TV j +oo. 
Proposition 2.4 shows that for e > satisfying (2.34) and (2.35), 7V(P-2)/2E|y(t, e)| 
uniformly in t G [TcT]. Proposition 2.5 proves that N(P-^y^E\V{t,e)Z-'^{t)\ uniformly 
in [To,T] for all e > 0. Then 

lim E[ sup \N^P-^^/^Fb(< MN{t) > -tN^^-P^/^EeMiA^ \] = 0- 

ATT+oo T„<t<T "-'^ 

Then supo<t<T |A^^^"^^/^i^b(< MAr(t) > -tN^^-P^/^K£,P)\ converges to zero in probability, 
since the probability of the events A^j^ can be made arbitrarily close to 1 by (2.21). This 
implies (2.33) and the proof of the lemma is complete. 

Proposition 2. 4: Assume that T > satisfies (2.10). Let us fix Q < e < 1/2 such that 
(2.34) and (2.35) hold. Then for any Tq > 0, Tq < T: 

lim Ar(P-2)/2E|Viv(«,e)| = (2.36) 

JVf+oo 
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uniformly in t E [To,T]. 

Proof. Let us estimate N^P~'^^/^E\VN{t, e)\ by the Cauchy-Schwartz inequality as in the 
proof of Proposition 2.2 for iV(P~^)/^E|T//v(i)|- After that we split it into four terms: 

iV(P-2)/2E|F^(t,e)| <[EWN{t,e)]'/^ = [lj,{t,e) - ll{t,e) +I%{t,e) + I%{t,e)Y/\ 

where 

WN{t,e) =m.,,,,y,(N^/^u(Ri,{a,a')) -Ei") (iV^^Vp (i?iv(cT, a")) -W 

,Hn {t,(T)+HN {t,CT')+HN {t,a")-3tN/2 



X e 
X 1 



{HNit,a}<Ntil+e),HN(t,a')<Nt(l+e),HNit,(r")<Nt(l+€)} 



•nP+mP + mg<e2/4 



X F{a ■ a' = miN, a ■ a" = m2N, a' ■ a" = mgiV} 

X ^eHN{t,a)+HN{t,<7')+HN{t,cr")-3tN/2 



X P{cr • cr' = miAT, cr • cr" = m2N, a' ■ a" = mgiV} 

X E[e-'^'^ (t,o-)+i/iv (t,o-')+-H'Ar {t,a")-3Nt/2 

^ I{-H'jv(t,T)>JVt(l + e) or /i'A,(t,cr')>Art(l + e), or H n {t,cr")>Nt{l+e)}] 

I%{t, e)= E (iV/j,(mi) - iv(2-p)/2E^f) (iV/p(m2) - iV^^-p)/^^^) 

e2/4<m^+mP+mg<3e 

X P{c7 • a' = miN, a ■ a" = m2N, a' ■ a" = maN} 

X ^^e^N{t,<7)+HN{t,a')+HN{t,a")-3Nt/2 

X ^{HN{t,CT)<Nt(l+e),HNit,<T')<m{l+e),HN{t,CT")<Nt{l+e)}] 

I%it,e)= Yl {Nfp{m,)-N^P-^y^Ee){Nfpim2)-N^P-^y^Ee) 



X P{(7 • cr' = miiV, a • a" = m2N, a' ■ a" = rriaN} 

X -^^e^N{t,a)+HN{t,a')+HN{t,a")-3tN/2 
X T^{HN(t,a)>m{l+e) OV HN{t,a')>Nt{l+e) Or fl"jv (t,CT")> Aft(H-e)}] • 

We will prove the uniform convergence to zero in [Tq, T] as iV f +oo of all these four terms. 
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The first term ljf(t) is not truncated and it refers to the configurations of spins with small 
correlations nii, m2 and 777.3. The proof of its uniform convergence to zero in [TcT] relies 
on (2.34) and it is completely analogous to the proof of the uniform convergence to zero of 
the sum l]^ + I%{t) + I%{t) in the proof of Proposition 1. Therefore, we omit the details. 

The second term I%{t) also contains only configurations of spins with very small cor- 
relations. If these correlations were zero, i. e. if HN{t,a) HN{t,a') and HN{t,a") were 
independent, then, indeed, the expectation involved in this term satisfies 

j^^^H^{t,.)+H^(ty)+H^it,a")-3Nt/2j^^^^ < 3E[e^«-^*/2 1^^^^^^^^^^] < exp{-iVieV2} 

is a standard Gaussian) by a well-known estimate for Gaussian random variables (5.1). 
We show that very small correlations 777 1, 7772, 7773 do not destroy the exponential convergence 
to zero of the corresponding expectation. Considering the third term I%{t), we neglect the 
truncation and use the asymptotic expansion (2.27) and condition (2.34). So we prove that 
the expectation ]Ee-^'v(t,o-)+i?iv(t,o-')+i?iv(t,o-")-3Art/2 multiplied by the probability of any given 

correlations goes to zero exponentially fast. Finally I%{t) refers to the configurations of 
spins with rather big correlations. Here, applying the estimate (5.1), we benefit from the 
truncation. The choice of e > according to (2.35) plays a crucial role in the analysis of this 
term. (Remember that this choice was possible only for T satisfying (2.10)). 

Now we proceed with the detailed proof. To treat the second term I%{t, e), we write 

' ii-{HN{t,a)>m{l+e) or HN{t,a')>Nt{l+e) or HN{t,a")>Nt{l+e)}\ 

_ 1pr ViVt(6+«2+€3)-3JVt/2 -(T -1 

— '^[^ ■ll-{gj>VlVt(l + e)Or ^2>\/lVt(H-e) or e3>V^(l + e)}J' 

where ,^1, ^2 and ^3 arc Gaussian random variables with zero mean, variance 1 and covari- 
ances cov(^i,^2) = /p(-Rjv(o-, o"')) = "^ii cov(^i,^3) = /p(i?jv(a, a")) = 777^, cov(^2,C3) = 
fp{RN{(j' ,a")) = ml, 777? + 777^ + 777^ < eV4. One gets 

=e^*.-3ivV2E[eVlvI(i+.)e.2^^^^^^^^^^g, 

where 7 = 1 + 7773 — (777? + 7773)^/2, /7 = 777? + 7773. Since 777? + 7773 < e^/4 < e, we may use 
the estimate for standard Gaussian random variables (5.1). It implies 

E[eV^(€i+«^+«3)-3ivt/22^^^^^^^^^^^] <C, exp{iVt(77^f + + 777^ - (e - 777? + 777^)V2)} 

<Ci exp{-ArroeV8} 
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for some constant Ci > 0, all t G [Tq, T] and all iV > 0, if + + < e'^/A, < e < 1/2. 
Thus 

sup E[e^^ (t,o-)+iJjv (t,o-)+iJjv {t,cr")-3Nt/2 

0<m^'+mf+mP<e2/4 

^ I{fl'iv(t,o-)>Art(H-e) or fl'iv(t,o-')>iVt(H-e) Or HN{t,(T")>Nt{l+e)}] 

< 3Ci exp{-iVroeV8} 

for all t E [To,T]. Since the other terms in Ij^it^e) have polynomial growth, the uniform 
convergence /^(t, e) — in [TqjT] follows. 

Let us turn to lf^{t,e). By the expansion (2.27) and condition (2.35) 

sup ^^^HN{t,a)+HN{t,a')+HN{t,a")-3Nt/2 
e2/4<m^'+mP+mf<3e 

X ^{HN{t,(T)<Nt(l+e),HN(t,(T')<Nt(l+e),HN{t,a")<Nt(l+e)}] 

X P{c7 • a' = miN, a ■ a" = m2N, a ■ a" = rusN} 

< C2 sup exp{N[t{m^ + + mf) — /(m)]} 

e2/4<m^'+mP+mf<3e 

< C2 sup exp{-/niV(m? + + m^)} < C2 exp{-/ne^/PiV/4} 

mP+mP+mg>e2/4 

for all t G [Tb,r], where C2 > 0, /ii > are constants. All other terms in 7|f(t, e) have 
polynomial growth, hence I%{t, e) — uniformly in [Tq, T]. 

Finally, consider lff{t,e). We have 

Iff \^HN{t.a)+HN{t,<T')+HN{t,(T")-3Nt/2 -ir 

ii^ e ' ii.{HN{t,a)<Nt{l+e),HN{t,a')<m{l+e),HN(t,a")<{l+e)Nt} 



< E 



^■^y Nt(3+2ml+2m^+2m^) i-3Nt/2j 



where ^ is a standard Gaussian, mi = /p(i?jv(o") (^')), 1^2 = /p(-RAr(c, (^")), 1^3 = /p(-Rjv(c', cr")). 
Since + + > 3e, we may apply the estimate (5.2). It yields 



E 



^HN{t,a)+HN{t,a')+HN{t,a")-3Nt/2 -ir 

< C3 exp{NtQp{ini,m2, ma, e)}. 



for some constant C3 > 0, alH G [Tq,T], N > and mf + m| + > 3e. On the other hand, 
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we also have: 

''^[^ '^{Hr,-{t,a)<Nt{l+e),HN{t,f7')<Nt{l+e),HN{t,<7")<{l+e)Nt} 
<r W L-H'jv(t,<T)+H;v(t,<T')+-H"jv(t,<T")-3Art/2 IT 

<^[e™^'^ ]i{i/^(t,<7')<Art(l+e),//jv(t,<T")<Aft(l+e)} 
_ F[p^«2+\/M53-3Art/2Tp/ V^€i I -J w 

-^[^ K2'«jJl{^2<VJVt(l+e),6<^/iVt(l+e)} 
gViVt(H-/:^2)6+^/^(l+M3)C3 j 



< e 



-Nt+Nta 



E 



{S2<VNt{l+e),^3<Vm{l+e)} 
g^JVt((l+M2)2 + (l+M3)"+2mP(l+M2)(l+M3))« J 



{-^2+2mg ^<2ViVt(H-e)} 



where , ^2 , ^3 are the same as in the analysis of the second term, ^ is standard Gaussian 
and 

a =(2m^ 7712^3 — m^^ — m^)/{2 — 2rni^) 
fi2 =(?TT'f — m2m^)/{l — mg^) 
Ms =("i2 - m^m|)/(l - mg^). 

One checks that 

+ /X2)2 + (1 + /X3)2 + 2mf (1 + /X2)(l + /xa) 

^ 2(1 + mg + (mf + m^)/2) ^ 2(1 + 3e/2) 
72 + 2mf ~ 72 + 2mf ' 

when + + > 3e. So, we are again in the position to apply (5.2). This yields 

W \ ^HN{t,a) + HN(t,a') + HN{t,cr")-SNt/2 -it 

^ ^ ' ^ ^HN{t,a)<Nt{l+e),HN{t,a')<Nt{l+e),HN{t,a")<{l+e)Nt} 

< C4 exp{tNLp{mi,m2, ma, e)}, 

where C4 > is a constant. Permuting mi , m2 and ma , we can derive in the same way that the 
same expectation does not exceed exp{tiVLp(mi, ma, m2, e)} and exp{fArLp(m2, ma, mi, e)} 
multiplied by some constant. Thus, taking into account (2.32), we obtain 

sup ^^^Hr^{t,a)+HN{t,a')+HN(t,a")-3Nt/2 
mj+m2+m3>ae 

X '^{HN(t,<7)<Nt{l+e),HN(t,cr')<Nt{l+e),HN(t,cr")<Nt{l+e)}] 

X P{cr • a' = miN, a ■ a" = m^N, a ■ a" = msN} 
< sup C5\/iVexp-! tArmin[(5p(mi,m2,ma, e), Lp(mi,m2,m3, e), 

mj+m2+m3>ae 

Lp(mi,ma,m2,e), Lp(m2, ma, mi, e)] - Ar/(mi, m2, m2)| 

(2.37) 

for all t G [0,ro], where C5 > is a constant. Now the relevance of the assumption (2.35) 
becomes clear. Due to (2.35), the right-hand side of (2.37) tends to zero exponentially fast. 
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as one can estimate it by C^V^ exp{— /i2iV}. The other terms in lfj{t,e) have polynomial 
growth. Thus I%{t, e) j +00 uniformly in [Tq, T\. This concludes the proof of the proposition. 

Proposition 2.5: For allT > satisfying (2.10) and all e> Q 

lim iV(P-2)/2E|V'jv(i,e)Zj^2(t)| =0 (2.38) 

ATf+oo 

uniformly in any interval [To,T], where <Tq <T. 
Proof. It follows from the definition of Vn {t) that 

N^P-^y^E\VN{t,e)Z-\t)\ < CNE^ J^HltT'''^'^'^^ (2-39) 

for all t > 0, where > is a constant. We will show that the expectation of this last 
fraction tends to zero exponentially fast. First of all, we observe that by (5.1) 

^^^^H^{t,a) - MH^{t,a)>m{l+e)} <& ' ■ (2-40) 

Let us represent the fraction in the right-hand side of (2.39) as 
Ifr pHN{t,(T) 

" ('2 411 

exp{lnE^e^^(*'<^) - ElnE<,e^'v(t,<T) + ]EinE<,e^'v(t,<T) _iVt/2}' 
To expand this formula, we will use the concentration of measure as in (5.3). The random 
variable Eae^^'^^^"'^ has the same distribution as 0(Ji, . . . , Jjvp), where the function 

^(xi,... jXatp) = InE^expjVtATi-P ^ Xi^M,- ,iv^iv^ii- ' ' 

is defined on Z'^'', Ji,... , Jatp are standard Gaussian random variables. The Lipschitz 
constant of 0(xi, ... ,Xf^p) is at most VtN^^VNP = \/tN . Substituting this function and 
u = Nte'^/4 into (5.3), we derive: 

P{| InE^.e^'^^*'^) - ElnE<,e-^^(*'^)| > Nte^/A} < exp{-Arte''/32}. (2.42) 

Let us introduce the events O^^ := {| In E<^e^^ - E In E^-e-^^^ > Nte^/4}. Conse- 
quently by (2.41) and (2.42) 

E^e ^-(''^H{H.(^t,.)>Ntil+e)} 

I,o^jE.e^-(*--)-^V^I,^,,,.)>^,,,,, (2 43) 

exp{lnE<,e-f^^(*'^) -ElnE^e-f^^(*''^) +ElnE<,e-f^'v(t,a) _^^/2} 

exp{ElnE<,e^^(*'<^) - Ari/2} 



E- 
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Observe that for any T satisfying (2.10) and any < Tq < T, there exists a constant > 
such that 

-K^/N <ElnE^e^''^^''^^ -Nt/2<0 (2.44) 

for all t G [IqjT]. The upper bound in (2.44) is immediate by Jensen inequality. Whenever 
the second moment of Ziq{t) truncated is finite, the left-hand side of (2.44) was established 
by Talagrand [Tl] in the analysis of the critical temperature. We will outline his proof in 
our situation. For given T satisfying (2.10), let us fix ?> 0, such that (2.34) and (2.35) hold. 
Let us define after 

By (5.2) there exists a constant i^i < such that 

EZjv(t,e) > i^i (2.45) 
for all i G [To,T]. Moreover, there exists a constant K2 > such that 

EZ%{t,Z)<K2 (2.46) 

for all t G [TcT]. The proof of (2.46) is analogous to the proof of the uniform convergence 
to zero of W"jv(i, e) in Proposition 2. We decompose Zis[{t,e) into four terms like it was for 
W"jv(i, e)- The last three of them go to zero uniformly in t G [To,r] and exponentially fast 
by the same arguments as I%{t), I%{t) and I%{t) do. We work out the first term similarly 
to the sum + 12 {t) + {t) in Proposition 1. The only difference is that tends to the 
integral along of the density of three independent standard Gaussians, which equals 1. 
Thus, in fact, ZN{t,'e) converges to 1 uniformly in [TcT] and (2.46) is obvious. Hence, for 
all t G [To,T] 

^^MM)^<M|M)!^<C:=;f, (2.47) 



(EZ„(t,£l)' (EZ„(t,£))' Kl 

Then starting from the Paley-Zygmund inequality and finally applying the concentration of 
measure inequality (5.3) with u = Nt/2 - ElnE^e''^^^*'^) + ln(ifi/2), we deduce 

I/4K3 < ' ' „ < ¥{ZNit,Z) > EZiv(i,e)/2} < P{E<,e'^^(*''^) > Kie^'/^/2} 

4{EZN{t,7)) 

= P{lnE^e^'^^*''^) - ElnE^e^~(*'^) > Nt/2 - E In E^e^^ + ln(iv:i/2)} 
< exp{[Ari/2-ElnE^e^^(*''") + ln(Ki/2)]V2A^O> 
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from where (2.44) follows. Finally, (2.40), (2.43) and (2.44) together imply 

and the proposition is proved.^ 

Proof of (1.19). To complete the proof of Theorem 1.3, it remains to show that 

lim inf Fp(mi,m2,m3) = 21n2. (2.49) 
After elaborating the functions Sp{m-\_,m2,ms) and Rp{m-\_,m2,ms), we get: 



Up(mi,m2,m3) =/(mi, m2, rn3)(l + mg) 

"(4(l + mg+ : ] + 



2 J + (1 - ) J (2.50) 

2(1 -mf) ^ 2 i^ + ^aJJ • 

It follows from (2.50) that for any p = 2k > 2 and any sequence {rrii^n, fn2,nj'n^i,n) £ -A. such 
that mi^n 1) 1^2,71 1, 'TT-a.n 1, as ri oo, 

lim yp(mi,„,m2,„,m3,„) = 21n2. (2.51) 

nt+cx> 

(In fact, by the definition of A we have ||mi| — |m2|| < 1 — l^nal for all {mi, 1712,1^3) G A, 
whence (m?_„ - „)2 = o(l - „).) Thus 

lim sup inf l^(mi, 777,2, ms) < 2 In 2. 
This fact and the next Proposition 2.6 together imply (2.49). <0> 

Proposition 2.6:Let {Pn} be a sequence of positive even numbers, Pn T +00. Assume that 
the sequence (mi_„, m2,n) '^3,n) € A satisfies one of the following conditions: 

(^nj i/iere e^si 5 > and a pair i and j, i,j = 1,2,3, i 7^ j, such that |mj^„| — ^ 1 and 
|^j,n| < 1 — (5 for all sufficiently large n; 
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(in) there exists S > such that |mi^„| <1 — S, \m2,n\ < 1 — ^; |?Ti3,„| < 1 — ^ for all sufficiently 
large n. Then 

lim inf yp„(mi,„,m2,„,m3,„) > 21n2. (2.52) 

nt+oo 

Proof: In the cases (i) and (iii) it suffices to substitute the sequence (w-i.n, 'm.2.n) ^3,n) into 
the function /(mi, m2, ma) (2/3 + (m^ + m2 +"1,3)""'^ ). Incase (ii) assume that e. g. |m3.„| 1 
and |mi^n| < 1 — <5. Then mf"^ = o(l). By definition of the set A we obtain | |mi^„| — |m2,n|| < 
1 — I"i3,n| ^ as n t +(X), thus = o(l) and (?77,f"„ — m2"„)^/(l — mg"^) = o(l). Moreover, 
if fnz,n 1) then mi^„ — m2,n and if m3^„ — > —1, then mi^„ + m2,n and therefore 
in both of these cases liminfjj^+00 -f("ii,n, m2,n, "i3,n) > ln2. This yields 
hm inf l^p„ (mi,^, m2,„, m3,„) 

nt+oo 

1 + ml\ 

> hm inf i7p„ (mi,„, m2,„, m3,„) > hm inf In 2 ' ' > 21n2 

nt+oo nt+oo ^ + 0\1.) 

and the proposition is proved. <C> 

3. The fluctuations of the partition function in the REM. 

Amazingly enough, the simplest of all our models, the REM, will be seen to offer in some 
sense the most interesting behaviour with regard to the fluctuations of the free energy. The 
main surprise here will be the existence of an intermediate region of temperatures where a 
CLT does not hold, but there a non-standard limit theorem will be proven. 

We begin with the proof of (i) of Theorem 1.4. 

Proposition 3.1: Whenever < /3 < ^JW1|2, 



Proof. This result will follow from the standard CLT for triangular arrays. Let us first write 

In = Inf 1 + . (3.2) 

We will show that the second term in the logarithm properly normalized will converge to a 
normal random variable. To see this, write 

Z^,;v-EZff,^r ^ ^ g-iv(i„2+/3V2) _giV/3V2^ ^ ^ y^^^y (3.3) 
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Note that E3^iv(a) = and ^yj^ia) = e-^(2in2-/32)[i _ g-iv/s^j ^nd thus 

Therefore we can write 

= .-f (.-.'.yiTF^^-l, 5: (3.5) 

where 3^Ar((T) = e^^'^^^'^~^^^[l — e~^^^]~^/^^iv(c) has mean zero and variance one. By the 
CLT for triangular arrays (see [Shi]), it follows readily that 

^ E yN{a)^M{0,l) (3.6) 
if the Lindeberg condition holds, that is in this case if for any e > 0, 

But 

oo 



1 



CX) 

, / e~^ci2; + o(l). 

(3.8) 

It is easy to check that the latter integral converges to zero if and only if /J^ < In 2/2. Using 
now the fact that = l+x + o{x) as a; — 0, it is now a trivial matter to deduce the assertion 
of the proposition. <0> 

Since the Lindeberg condition clearly fails for 2/3^ > In 2, it is clear that we cannot expect a 
simple CLT beyond this regime. Such a failure of a CLT is always a problem related to "heavy 
tails", and results from the fact that extremal events begin to influence the fluctuations of 
the sum. It appears therefore reasonable to separate form the sum the terms where Xa- is 
anomalously large. For Gaussian r.v.'s it is well known that the right scale of separation is 
given by un{x) defined by 



oo 



>N I ^e-.V2 = e-x (39) 



/27r 

un{x) 
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which (for x > — lnA^/ln2) is equal to (see e.g. [LLR]) 

x G M is a parameter. Let us now define 

Z%^^ ^ E.e^^^<'I{x<,<„.(.)}. (3.11) 

We may write 

Zf3,N '^0,N ~ ^^P,N _|_ ^I3,N - -^g^AT - ^(Zff^jv - -^g^jy) 

Let us first consider the last summand. We introduce the random variable 

y^Nix) = ^^1^^ = e--0n^+^V2) J2 e'^^--^I{x.>..(.)> (3.13) 

It will be convenient to rewrite this as (we ignore the subleading corrections to un{x) and 
only keep the explicit representation (3.10)) 



>V.(x)=e--(-^+^V2) J2 e/^^-(.-(x.))j^^_^^^^^^^^^ 



Let us now introduce the point process on M given by 

A classical result from the theory of extreme order statistics (sec e.g. [LLR]) asserts that the 
point process Vn converges weakly to a Poisson point process on M with intensity measure 
e~^dx. We can, of course, write 



E ^^""'^'''^%-(x.)>.} = / e-^VM{dz), (3.16) 



where we set a = /3/\/21n2. Clearly, the weak convergence of Vn to V implies convergence 
in law of the right hand side of (3.16), provided that e"^ is integrable on [x, oo) w.r.t. the 
Poisson process with intensity e~^. This is, in fact never a problem: the Poisson point 
process has almost surely support on a finite set, and therefore e"^ always a.s. integrable. 
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Note, however, that for P > \/2 In 2 the mean of the integral is infinite, indicating the passage 
to the low temperature regime. Note also that the variance of the integral is finite exactly 
if a < 1/2, i.e. < In 2/2, i.e. when the CLT holds. On the other hand, the mean of the 
integral diverges if x J, oo; note that at minus infinity the points of the Poisson point process 
accumulate, and there is no finite support argument as before that would assure the existence 
if X is taken to — oo. The following lemma provides the first step in the proof of part (ii) of 
Theorem 1.4 and of Theorem 1.5: 

Lemma 3.2: Let yVN{x),a be defined above, and let V be the Poisson point process with 
intensity measure e~^dz. Then 

CO 

gf(y2hr2-/3)^ + t[ln(iVln2)+ln4.]-^^(^)^ j^<^z^i^dz). (3.17) 

X 



Remark: Note that the mean of the right hand side is finite if and only of /3 < \/21n2. Thus 
only in that case does this lemma also allow to deal with the centered variable appearing in 
(3.12). 



We now need to turn to the remaining term. 



I3,N 



Vn{x) 



(3.18) 
(3.19) 



where 

Viv(x)^Z|_jv-EZ|,iv 

One might first hope that this term upon proper scaling would converge to a Gaussian; 
however, one can easily check that this is not the case (the Lindeberg condition will not be 
verified). However, it will not be hard to compute all moments of this term: 

liemma. 3.3: Let Vn{x) be defined by (3.19). Then for a > 1/2 and any integer k > 2 



lim 



2- ArgJV/3v/2hr2- f [ln(Ar In 2)+ln 47r] 



^(ka — i)x 



i\ ^ tx\...ii\{i^a-\)...{lia-\) 



For a = 1/2, we have for k even 



lim 

ATT+oo 



E[VNix)f 



A;! 



2-AfgAf/3V21n2 



_ (fc- 1)!! 

Jk/2)\¥ ~ 2^/2 



(3.20) 



(3.21) 
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and for k odd 

JVT+oo ^ ^ 



2-JVgJV/9V21n2 

(which are the moments of the normal distribution with variance 1/2). 

Proof. This is a pure computation. Set rjv(cr) = e^^^" I{x<,<ujv(x)}- Note that for 

ET^r(a)= / ^e-^+^^^ = e^'3'/2fl- / ^e'^) ~ e'^'^/^ (3.23) 
7 \/27r ^ J yl'K ' 

while for /? > V21n2 and ah > 1, and for /? > A/ln2/2 and for > 2, 

ttAr(a;) Wat (a;) — fe/3-\/JV 

-i -i (3.24) 

^ gJVfc2/3^/2_£ ^ f__^_gfe[/3^/2h^2A^+Q;x-f [ln(Arin2)+ln47r]] 

V2^{kpVN - un{x)) ka-1 
Formula (3.24) is also valid for /3 = v'21n2 with A; > 1 and for /3 = A/ln2/2 with A; > 2. 
It is easy to see from the computations above that for /3 = \/21n2 with A; = 1 and also for 
P = y/\n2/2 with A; = 2 we have 

E[T^(C7)]'= - = efc[/3V2hr2iV+ax]_ (3 35) 



We set r7v(cr) = 2-^riv(cr); by (3.24) we get for (3 > ^\a.2/2 with A; > 2 and also for 
(3 > V21n2 with k>l 



E[Tiv((T)]'' = ^ ^ ^fc[/3V2hr2Ar-ln2+ax-^fln(JVln2)+ln47r]]_ (■g^g) 



This formula is also true for /3 = A/ln2/2, /c > 2 and /3 = V21n2, A; > 1. For /3 = V2 In 2 
and A; = 1 and also for j3 = Y^ln2/2 and fc = 2 by (3.25) 

Now 

E [VN{x)f = e( ^ [rAr(a) - ETiv(c7)]J = ^ E J][rjv(ai) - ETiv(ai)] 

ct&Sn cri,...,ah&SN »=1 

i=l <i,...,<!j>2 1 »■ V / 

(3.28) 
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Note finally that for Z > 2 and /3 > ^J\n.2/2 



Section 3 



E[rjv(a) - ET^rCa)]^ = Y^-^)' QErjv(c7)^-^ [Erjv(c7)p ~ ETjv(a)^ (3.29) 

In fact, if ^^n272 < ^ < ^21^, Z > 2, j > 1, j 7^ Z - 1, Z, then by (3.23) and (3.26), (3.27) 

E[r^-^(^)][ETjv(a)p' Nj{pV2-pV2T^)^o ( AT^i/s^ .3 oqn 

E[T^(a)] ^ V ; ^ • ^ 

For Z>2, j = Z-l,Z 

E[rj,(a)] V ; - ^ 



For P > V2h^, Z > 2 and i > 1 by (3.26) and (3.27) 



E[T^>)l[Er.(.)p^ (3.32) 
Thus for Z > 2 and /? > A/ln2/2 and also for Z > 3 and /3 = ^^\n2/2 

r) — N —X 

ErTW(a) - ETW(a)]^ = , ^ r2-iVgJV/3V2hr2gaXg-f [ln(iVl„2)+ln4.]1^_ .3 33) 

A;a — 1 

Inserting this result into (3.28) gives the assertion of the lemma (3.20). 
For (3 = v^ln2/2 and Z = 2 by (3.27) we have 

Q — TV — X 

E[TW(a) -ETW(a)]' = ^ (3.34) 

Inserting this formula into (3.28) we see, that the term with li, . . . ,1^ = 2, i = k/2 brings the 
main contribution to the sum, and all others are of smaller order, because of the polynomial 
terms g-'f '"(^'"2) (3 33) This implies (3.21) and (3.22) and the lemma is proved. 

Remark: One sees that if we let x | —00, and rescale properly, the corresponding moments 
converge to that of a centered Gaussian r.v. This could alternatively be seen by checking 
that the Lindeberg condition holds for the truncated variables provided x < — 21nln2^. 

A standard consequence of Lemma 3.3 is the weak convergence of the normalized version 
ofVjv(x): 
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Corollary 3.4: For ^\n2/2 < (i, 

(^2hr2-/3)^ + tMiVln2)+ln4.] Viv(^ V ^^^^^^^ ^3 3^^ 

where V(a;, a) is the random variable with mean zero and kth moments given by the right 
hand side of (3.20). For (3 = ^\n2/2 

^Ar(0,l). (3.36) 

The next proposition will imply (ii) of Theorem 1.4. 

Proposition 3.5: Let y^ln2/2 < f3 < \/2 In 2. Then for a; G M chosen arbitrarily, 

oo oo 

(V2hr2-/3)^+t[i„(ivi„2)+i„4.]i^^^ ^ y^^^^^ ^ /e°^P(dz) - / e«^e-^dz, (3.37) 

X X 

where V{x,a) and V are independent random variables. 

Proof. (3.37) would be immediate from Lemma 3.2 and Corollary 3.4, if WAr(a;) and Vn{x) 
were independent. However, while this is not true, they are not far from independent. To 
see this, note that if we condition on the number of variables X^, njv(a;), that exceed un{x), 
the decomposition in (3.12) is independent. On the other hand, one readily verifies that 
Corollary 3.4 also holds under the conditional law P[-|n7v(a;) = n], for any finite n, with the 
same right hand side V{x,a). But this implies that the limit can be written as the sum of 
two independent random variables, as desired. <C> 

Since for /J^ > ln2/2, a > 1/2, one sees that EV{x,af = e^^^^-^VC^a - 1) tends to zero 
as a; 1 — oo. Therefore we see that 

X X 

V{x,a)=T:) lim [ e'^^Pidz) - [ e'^'e'^dz (3.38) 

2/T+oo J J 

-y -y 

which means that wc can give sense to the Poisson integral e°'^{V{dz) — e~^dz) We see 
that Propositions 3.1 and 3.5 imply Theorem 1.4. 

Remark: The appearance of the intermediate region with non-Gaussian fluctuations may 
appear surprising in view of the fact that in the p-spin models, we could prove the CLT up 
to a much higher value of /3, in fact up to almost the critical value. The reason, however, 
lies in the fact that in the p-spin model the Gaussian part of the fluctuation is always on a 
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polynomial scale in N, while the truncation error ((Z^^^r — ZjjY)/EZ^^jv) is exponentially 
small even when we truncate at + e)\/iV, way below where we truncate in the REM. This 
means that the CLT contribution will always dominate the extremal fluctuations. In the 
REM everything is exponentially small, and while a sufficiently truncated partition function 
gives a Gaussian contribution, this is dominated by the larger extremal fluctuations in the 
intermediate regime. In other words, the extra correlations in the p-spin models strengthen 
the Gaussian fluctuations more than the extremal ones which sounds intuitive. 

We now turn to the 

Proof of Theorem 1.5. We will see that the computions above almost suffice to conclude 
the low temperature case as well. With the notations from above, we write 

Z(3,N = -^|,jv + {^P,N - Z^^n) (3-39) 

Clearly for (5 > V21n2 

Zp^I, - = e^[/3V2hr2-I„2]-f [ln(iVl„2)+l„4.] ^ ^^^^^^^au-\X^) ^g^g) 

SO that for any a; G M, 

oo 

{Zp^N - Z-^)e-^[/3^2h^-l°2] + f [ln(iVl„2)+l„4.] ^ j ^^^-p^dz). (3.41) 

X 

Now write 



ZIn = ^ZIn i 1 + ^''^^f^"'"' ) • (3-42) 



Let us first treat the case /? > \/21n2. By (3.24) we have 



EZ^ ~ f__^_g/3V2hr2Ar+aa;-f [ln(iVln2)+ln47r] ^3 ^g-j 

a — 1 



Thus 

7x 



-iV[/3V2hr2-ln2] + f[ln(iVl„2)+ln4.]^x^ ^ + ^IliL-^:!!:^) (1 + o(l)). (3.44) 



IE2' x((x — i) 

Using Lemma 3.3 we see that now '^'Z^-^ — — _-. converges in distribution to a random 
variable with moments given by the right hand side of (3.20). Moreover, as x J, —00, this 
variable converges to zero in probability. Since the same is true for the pref actor, the assertion 
of the theorem is now immediate. 
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Let US consider now the case P = \/21n2. Proceeding as in (3.24), 



«Ar(0)-V2JV In 2 

— CX) 

We use the decomposition 

Z0,N = Zn^N - Z^^,N + ^Z^^j^ + {Z'^^j^ - lEZ^^jv). (3.46) 
By (3.45), EZ^ j^/EZp^N ~ 1/2. By (3.14), we see easily that 

^M_3^ = Wn{x)^0 a.s. (3.47) 
even though E>VAr(0) = 1/2! Thus the more precise statement consists in saying that 

oo 


Note that of course the hmiting varaible has infinite mean, but is a.s. finite. Finally, by 
Corollary 3.4, 

EZis^N 

The same arguments as those given after Proposition 3.5 allow us to identify V(0, 1) with 
the Centered Poisson integral i'Pidz) — e~^dz) . This implies (1.24). (1.25) is an 

immediate corollary. This concludes the proof of Theorem 1.5.<0><O 

Appendix 1. Some remarks on the case p odd 

Conjecture 4.1:Let p = 2k + 1, k > 1. There exists (3p > such that for all (3 < Pp 

^(p-2)/2i^^^^^(^) (4.1) 

in distribution as N 1 +oo, where Moo{t) is a centred Gaussian process with independent 
increments and 

Moreover f3p — V2 In 2, as p\ +oo . 

Discussion. Let us try to adapt the martingale method in this case. This leads to 



VN{t) =E„,„>[N{RM{a,a')Y - N^-HE^^) 



^HN{t,a)+HN{t,a')-Nt _ 
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Then 

NP-^EVN{t) = Yl [NP-^mP - tE^2p^ e*^'">(a • a' = mN). (4.2) 

m=0,±l/N,... ,±1 

It is easy to show that A^^~^EVAr(t) ^ as | +00 for all t such that t < info<m<i 4'{m)m~P . 
As in the proof for p even, we can concentrate only on configurations of spins with corre- 
lations ra close to zero, since others bring an exponentially small contribution. Note that 
P((j • a' = mN) = P((7 • a' = —mN) and consequently I{m) = I{-m) = — m^/2(l + o(l)), 
m — > 0. Summing up the terms in (4.2) with correlations m and — m, we get 



iVP-2Ey,v(i) =^2= NP-^mPie*^"^" - e-*^™^)e-^-f("') - 2tEfP + o(l) 

" 7T).>0, 

|m|<A'"-l/3-'5 

= ^ iVP-^mf (tiVmP) (1 + o{l))e-^^^"'^ - 2tE^2p ^ 

* ^ m>0, 

|m|<Arl/3-« 

00 

J s^Pe-'^/^ds - 2iE^2p + 0(1) ^0, AT t +00. 



27r 



Moreover, as for p even, it is also not difficult to show that the truncated value N^P~'^^'EVN{t, e) 
tends to zero for all t up to Talagrand's bound (2.16). 

Let us now try to perform a rigorous proof of Conjecture 4.1. Proceeding along the lines 
of the proof for p even, we come to the problem of convergence iVP~^E|Vjv(i)| 0. To get 
rid of the absolute value of VN{t), let us first apply the Cauchy-Schwartz inequality in the 
same way as it was in the proof of Proposition 2.2. We obtain 

[N^P-^^E\VN{t)\]^ < Yl i^^'^rrfi - tEe){NP-^m^2 - tE,e'')e^*^™^^"^^™^^ 

mi, 7712, ms (4-3) 

X P((7 • a' = miN, a-a" = m2N, a' ■ a" = mzN). 

Surprisingly, the right-hand side of (4.3) does not converge to zero. The problem arises from 
the fact that 

I{mx,m2,ms) = /(-mi, -m2, ma) = /(mi, -m2, -ms) = /(-mi, m2, -ma), 

but 

I{mi,m2,m3) / I{mi,m2, -ms). 

In fact, opening the brackets in {NP~^m\ — tE^'^P){NP~^m2 — tE^^^) one can split the right- 
hand side of (4.3) into four terms. Let us elaborate the first one summing up together 
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the terms with correlations, having the same absolute values \mi\, \m2\, {msl and the same 
probability: 

mi>0,m2>0,m3>0 
mi>0,m2>0,m3>0 
mi>0,m2>0,m3>0 

Jg-A'/(™i,m2,"i3)+JV(m^+m^+m^)/2 _ ^-N I (mi,m2,-m3) + N {mf+ml+ml) /2j 

+ t^{E^^Pf +o{l). 

This term is of order N'^p-^'^/H{E^p+^)^{1 + o(l)) + t^E^^p^ gj^^^g ^^le expansion of 

^^-NI{mi,m2,m3)+N{ml+m2+ml)/2 _ ^-NI{mi,m2,-m3)+N{ml+ml+ml)/2-^ 

the main term is of order Nmim2ms- The sum of all other three terms in (4.3) tends to 
-t2(E^2p)2 rpj^^g ^j^g right-hand side of (4.3) is of order N^-^y^tiES,P+'^)^ and it does not 
converge to zero for N | +oo. Therefore the proof for p even is not suitable at this point for 
p odd. 

A possible solution for this problem is to apply the Cauchy-Schwartz inequality in a dif- 
ferent way passing to the fourth moment of ^jv(i): 

[N^P-^^E\VN{t)\Y <EE,,a',a'',a"'(NP-'(^RN{a,a')) -tEe'')[^^^^) -tE^'') 

^ ^HN{t,<T)+HN{t,<7')+HN{t,a")+HN{t,<7"')-2Nt 

It can be proved that the right-hand side of this last inequality tends to zero for all t up to 
some bound. But technical details are very tedious. We will only say that six parameters 
mi, . . . ,mQ have to be considered. The group of 64 correlations with fixed absolute values 
|mi|, . . . , Imel splits into eight groups of correlations having the same probabilities. 

Furthermore, it will be technically even much harder to extend the bound of t by the 
truncation of the Hamiltonian. We will have to taJie into account five different cases and 
their permutations where some of correlations are large and some are small. Each of these 
cases will require very tough computations. 
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Appendix 2. Two useful theorems 

Proposition 5.1: Let ^ be a Gaussian random variable with = 0, E^^ = 1. Then for all 
a,b> 

E[exp{aOE{«>6}] < ^^^^^ _ exp{-b''/2 + ab}, ifb>a, (5.1) 

E[exp{aC}I{^<6}] < , j exp{-6V2 + ab}, ifb<a. (5.2) 

y/2Tr{a — b) 

Theorem 5.2: Assume that f{xi, . . . ,Xd) is a function on with a Lipschitz constant 
L. Let Ji, . . . ,Jd be independent standard Gaussian random variables. Then for any u > 

P{|/(Ji, ...,Jd)- E/(Ji, ...,Jd)\>u}< exp{-^.V(2^')}- (5-3) 
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